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Amplitude  and  Phase  Curves  for  Ground -Wave  Propagation 
in  the  Band  200  Cycles  Per  Second  to  500  Kilocycles 

_ 

James  R.  Wait  and  H.  Herbert  Howe  ^ 

After  making  several  extensions  to  the  formulas  of  Van  cler  Pol  and  Bremmer,  field 
strength  and  phase  values  of  the  very  low  frecpiency  ground  wave  from  a  short  vertical 
antenna  are  computed.  The  ground  conductivity  values  chosen  are  4,  0.01,  and  0.001  mho 
per  meter.  The  distances  considered  range  from  1  to  1,500  miles. 


1.  Introduction 

Considerable  interest  has  been  shown  recently 
in  the  propagation  characteristics  of  very  low 
frequency  radio  waves  [1,2,  3].^  Electromagnetic 
energy  is  present  in  lightning  discharges  at  fre¬ 
quencies  down  to  100  cps  and  less  [4].  There  has 
also  been  a  revived  interest  in  very  low  frequencies 
(VLF)  for  communication  purposes  because  they 
are  onh^  slightl}"  affected  by  ionospheric  disturb¬ 
ances  for  ranges  up  to  200  miles  or  so.  VLF  has 
been  found  particularly  suitable  for  navigational 
systems  using  a  phase-comparison  technique  [5,  6]. 
In  such  a  system,  two  continuous-wave  signals 
are  converted  in  the  receiver  to  a  common  com¬ 
parison  frequency,  which  is  displayed  on  a  phase 
meter.  Lines  of  equiphase-difference  form  the 
position  lines  from  which  this  navigational  aid, 
Imown  as  ^‘Decca”,  is  based.  Another  promising 
application  of  VLF  is  to  geophysical  prospecting, 
where  it  is  known  that  the  propagation  of  the 
ground  wave  is  dependent  on  the  upper  strata  of 
the  earth’s  crust.  At  frequencies  of  the  order  of 
100  kc,  the  substrata  at  depths  down  to  several 
hundred  feet  will  influence  the  attenuation  and 
phase  velocity  of  the  ground  wave  [7]. 

In  view  of  the  present  VLF  developments  in  the 
above-mentioned  fields,  it  seems  very  desirable  to 
present  detailed  information  on  the  ground-wave 
characteristics  of  an  electric  dipole  source.  The 
frequencies  considered  are  200  cps  to  500  kc. 
The  conductivity  between  the  source  and  the 
observer  is  considered  to  be  essentially  homogene¬ 
ous,  with  conductivity  c  in  mhos  per  meter.  The 
values  of  a  chosen  are  4,  10~^,  and  10~^,  cor¬ 
responding  to  sea  water,  well-conducting  land,  and 
poorly  conducting  land,  respectively.  In  this 
frequency  range  the  displacement  currents  in  the 
ground  can  be  neglected  compared  to  the  conduc¬ 
tion  current,  so  the  dielectric-constant  value  need 
not  be  considered.  The  validity  of  this  approxi¬ 
mation  was  examined  in  detail  elsewhere  [1,  8]. 


1  The  contribution  of  the  second  author  consisted  of  devising  the  program 
for  the  automatic  computer.  The  responsibility  for  the  technical  content  of 
this  paper  rests  completely  on  the  first  author. 

2  Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this  paper 


The  distances  in  miles  from  the  source  to  the 
receiver  are  taken  from  1  to  1,500,  inclusive.  To 
account  for  normal  atmospheric  refraction  [9], 
the  effective  earth  radius  is  taken  to  be  k  times  the 
actual  radius.  It  should  be  pointed  out,  however, 
that  the  effective  earth-radius  concept  is  not 
strictly  valid  at  very  low  radio  frequencies.  For 
most  of  the  calculations,  k  is  set  equal  to  4/3  to 
conform  with  standard  practice. 

2.  Some  Theoretical  Considerations 

For  any  usual  source  of  VLF  waves,  the  equiva¬ 
lent  antenna  will  have  a  height  that  is  small  com¬ 
pared  to  the  wavelength.  Furthermore,  for  most 
cases  the  distance  of  the  ol)server  to  the  source  will 
be  large  compared  to  the  height  of  antenna  or 
vertical  column  of  current.  In  this  case,  the 
fields  are  proportional  to  the  product  of  the  aver¬ 
age  current  /  and  the  height  I  of  the  antenna.  The 
solution  for  the  field  of  a  vertical  dipole  of  moment 
II  over  a  spherical  homogeneous  earth  was  first 
carried  out  successfully  by  Watson  [10].  Methods 
for  obtaijiing  numerical  values  for  application  to 
radio-wave  propagation  have  been  developed  by 
Van  der  Pol  and  Bremmer  [11],  Norton  [12],  and 
others.  The  application  of  these  methods  is  very 
straightforward  for  the  purpose  of  calculating  the 
amplitude  of  the  ground  wave  for  medium  and 
high  radiofrequencies.  Some  special  considera¬ 
tion  is  required,  however,  at  very  low  radiofre¬ 
quencies  where  the  induction  and  static  fields  of 
the  antenna  are  not  negligible.  This  is  particu¬ 
larly  true  from  the  standpoint  of  the  phase  varia¬ 
tions  with  distance. 

Choosing  a  spherical  coordinate  system  {r,d,4>), 
the  earth’s  surface  is  defined  by  r=a,  and  the 
source  is  a  radially  oriented  electric  dipole  of 
moment  II  located  on  the  earth  at  0  =  0.  A  time 
factor  expiioot)  is  imj)lied  throughout.  The  radial 
electric-held  component  Er  in  the  air  caii  be 
derived  from  a  scalar  function  n  by 
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The  vertical  component  E  of  the  electric  field 
at  a  distance  D  measured  along  the  surface  is  then 
given  by 


E=>- 


1  d 


D  dD 
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(2) 


since  e=D/a  is  small  compared  to  unity,  Van  der 
Pol  and  Bremmer  [11]  give  the  folio v\dng  expression 
lor  IT : 
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where  d=27r/free-space  wavelength,  and  eo= 
8.854X10“^^  farad/m.  TP  is  the  now  famous 
residue  series  given  by 


TP=  [27r(^a)i/3  1 /2g- ^  W4 
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and  where  the  Ts  are  solutions  of  the  equation 
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where  the  left-hand  side  is  the  ratio  of  two  Hankel 
functions  of  the  second  kind.  Van  der  Pol  ami 
Bremmer  develop  methods  for  obtaining  numerical 
values  of  r,  [3,  11]. 

The  ‘Residue  series”  representation  for  TP  is 
quite  suitable  for  computation  when  D  is  large 
enough  for  the  series  to  converge  rapidly.  At 
shorter  distances  (down  to  1  mile)  and  very  low 
radio  frequencies  the  series  would  become  im¬ 
practical,  since  thousands  of  terms  would  be  re¬ 
quired  to  secure  convergence.  At  these  verv 
short  distances,  however,  the  earth  can  be  con"- 
^i^d  another  expression  for 
H  {p)  IS  more  suitable,  which  is  given  by 


B  — 1  i(7rp)^2__2p+{(7r)'/2 


-j^(8-165'‘+7a'*)p3+ 


(6) 


yiere  p  is  .Sommerfeld’s  luimei-ieal  distance  do- 
hned  by  ’ 


f) 
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The  above  series  formula  for  TT"  was  proposed  by 
Bremmer.^  It  can  be  verified  that  when  6=0 
corresponding  to  a  flat  earth,  the  right-hand  side 
corresponds  to  the  series  expansion  for  the  well- 
known  Sommerfeld  flat-earth  propagation  func- 
tion  y(p),  which  has  now  been  tabulated  exten¬ 
sively  [12].  It  has  also  been  shown  recently  that 
the  general  formula  for  TP  can  be  extended  to  a 
^^^’‘"^dified  ground  with  elevated  antim  as 

Using  (2)  and  (3),  the  electric  field  is  now  given 
bv 
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where  (3  ^(1  + 1/2 13ak), 

The  three  principal  terms  in  the  above  expression 
proportional  to  l/D,  1//U,  and  1/D^  can  be  identi¬ 
fied  with  the  radiation,  induction,  and  static  fields 
respectively.  It  is  clear  from  the  above  develop¬ 
ment  that  the  effect  of  finite  conductivity  and 
earth  curvature  will  affect  to  some  extent  all  three 
components  of  the  field.  At  very  short  distances 
where  and  D/a  are  small  com¬ 

pared  to  one,  the  expression  for  the  field  simplifies 
considerabl}^  to 


g- 
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which,  of  course,  is  identical  to  the  field  of  an 
electric  dipole  on  a  flat  perfectly  conducting  plane. 
At  large  distances  where  (3D:^1,  the  expression 
lor  E  is  given  by 


2TrD 


(10) 


wliere  W  is  defined  by  eq  (4)  or  (6).  At  inter- 
meiliate  distances  whore  neitiier  of  the  above  two 
special  cases  is  valid,  it  appears  to  be  necessarv 
m  graerid,  to  employ  the  above  p  series  formula 
for  A.  For  most  purposes,  however,  it  is  suffi¬ 
ciently  accurate  to  use  the  following  formula  for 
the  whole  range  of  D. 


E= 


2TrD 


3  Private  communication. 
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The  maximum  error  in  this  formula  occurs  where 
is  in  the  vicinity  of  unity  and  for  the  range  of 
parameters  chosen  in  this  report  will  affect  the 
amplitude  to  within  1  percent  and  the  phase  by 
less  than  2  degrees.  Therefore,  eq  (11)  can  be 
used  for  all  the  calculations. 

3.  The  Numerical  Results 

The  major  numerical  task  in  this  work  was  the 
evaluation  of  the  residues  series  for  W.  In  appli¬ 
cations  to  field-strength  calculation  at  medium 
radiofrequencies  it  is  usually  sufficient  to  use  only 
one  term  of  the  residue  series  and  then  to  inter¬ 
polate  with  the  fiat-earth  formula,  which  is  valid  at 
short  distances.  The  procedure  for  making  rapid 
calculations  of  this  type  has  been  developed  b}^ 
Norton  [12].  For  amplitude  and  phase  calcula¬ 
tions  at  VLF,  the  interpolation  method  is  not  of 
sufficient  accuracy,  and  more  terms  of  the  residue 
series  are  required.  In  fact,  it  was  considered 


Table  1.  Amplitude  of  W  (k  =  4l3) 


Frequency 

Distance  (miles) 

37.7 

75.3 

150.6 

377 

753 

1,506 

(7  =  4 

kc 

0.2 

1.00000 

0.  99407 

0.  97581 

0.  93220 

0.  81890 

.5 

. 99900 

. 99047 

. 96176 

.  89458 

. 72915 

1.0 

.  99850 

. 98641 

.  94614 

. 85400 

.  64030 

2.0 

.  99779 

0. 99334 

.  98069 

.  92445 

.  79981 

. 53391 

5.0 

.  99638 

.  98931 

. 96939 

.  88286 

. 70267 

. 37563 

10.0 

.  99478 

. 98478 

. 95680 

.  83827 

.  60807 

. 25677 

20.0 

.  99252 

.  97839 

. 93925 

.  77923 

. 49714 

. 15426 

50.0 

.  98802 

. 96582 

. 90541 

.  67493 

. 33772 

. 06105 

100.0 

.  98296 

. 95183 

.  86886 

.  .57555 

.  22340 

. 02377 

200.0 

. 97578 

. 93237 

.  81996 

. 46216 

.  12951 

.00714 

500.0 

. 96130 

. 89459 

.73171 

. 30572 

.  04907 

. 00089 

<t  =  10“2 

0.2 

1.00060 

0.  99873 

0. 99401 

0.  97582 

0.  93223 

0.  81899 

.5 

. 99040 

.  96179 

.  89468 

. 72946 

1.0 

.  98634 

.  94620 

.  8542.5 

. 64104 

2.0 

’.'99843 

‘.’99333 

.  98062 

.92461 

. 80044 

. 53554 

5.0 

. 99914 

. 98927 

. 96930 

.88333 

.  70458 

.  379.52 

10.0 

. 99532 

. 95652 

. 83914 

.61196 

. 26304 

20.0 

.  99252 

. 97783 

.  93797 

. 77999 

.  .50364 

. 16225 

50.0 

.  98508 

. 95997 

. 89577 

. 66732 

.34198 

. 06646 

100.0 

.  96884 

. 92669 

. 82821 

. 53124 

. 20598 

. 02256 

200.  0 

.  92053 

.  83428 

. 66792 

. 30708 

. 06808 

. 00243 

500.  0 

— 

— 

— 

— 

. 00000 

<7  =  10-3 

0.2 

1.00060 

0.  99873 

0. 99434 

0.  97583 

0.  93228 

0.  81916 

.5 

. 99732 

. 99062 

. 96182 

.  89485 

.73010 

1.0 

.  99936 

.  99.572 

.  98650 

. 94625 

.  85465 

. 64246 

2.0 

. 99838 

.  99324 

.  980.56 

. 92458 

.80121 

. 53834 

5.0 

.  99653 

.  98866 

. 96837 

. 88212 

. 70525 

. 38440 

10.0 

. 99365 

. 98207 

. 95224 

. 83247 

.  60808 

. 26627 

20.0 

.  98687 

.  96739 

.92012 

. 75019 

. 47724 

. 15362 

50.0 

inn  (\ 

. 95220 

.  89869 

. 79018 

. 50525 

. 21268 

. 02962 

lUU.  U 

200.  0 

kkl’’ 

k  ’’’’ 

. 00091 

.  00000 

500.0 

.  02701 

. 00912 

’.’00069 

. 00001 

. 00000 

Distance  (kilometers) 

60.  6 

121 

242 

606 

1,211 

2,  420 

desirable  at  the  outset  of  this  work  to  carry  out 
calculations  from  the  residues  series  formula  at 
distances  down  to  37  miles,  where  no  interpolation 
with  the  p  series  formula  is  required.  In  fact, 
there  was  a  region  of  overlap  in  the  calculations 
that  provided  a  convenient  check  on  the  computa¬ 
tions.  Due  to  the  excessively  large  number  (up  to 
1,000)  of  terms  required  to  secure  convergence  of 
the  residue  series  at  sliort  distances  and  frequencies 
as  low  as  200  cps,  it  was  summed  on  the  SEAC  (the 
National  Bureau  of  Standards  electronic  auto¬ 
matic  computer).  Values  of  the  amplitude  and 
phase  of  11"  are  provided  in  the  tables  1  and  2, 
respectively,  for  several  conductivity  values  and 
various  distances.  There  are  several  gaps  in  the 
tables  that  correspond  to  situations  where  it  was 
not  possible  to  devise  a  convenient  program.  In 
some  cases  this  occurred  where  the  Bremmer-type 
[3]  series  formulas  for  5  did  not  converge. 

Table  2.  Phase  of  W  (k  —  4l3) 


(The  values  are  all  negative) 


Distance  (miles) 

F  requency 

37.7 

75.3 

150.6 

377 

753 

1,506 

<t  =  4 

kc 

0.2 

0.  050 

0.  364 

1.  431 

4.  021 

11.  Ill 

.5 

.079 

.579 

2.  267 

6.  323 

17.  222 

1.0 

.  112 

.  825 

3.209 

8.  883 

23.  829 

2.  0 

.  162 

0.  429 

1. 178 

4.  541 

12.  445 

.32.  717 

.5.0 

.273 

.  704 

1.898 

7. 183 

19.  .319 

49.  020 

10.  0 

.414 

1.036 

2.  738 

10.  160 

26.  792 

65.  744 

20.  0 

.641 

1.546 

3.  978 

14.  377 

36.  939 

87.  334 

50,  0 

1.  190 

2.  693 

6.  616 

22.  793 

55.  986 

126.  041 

100.  0 

1.965 

4.  205 

9.  879 

32.  423 

76.  410 

168.  527 

200.  0 

.3.  .341 

6.  740 

15.  031 

46.  466 

104.  711 

222. 115 

500.0 

7.  052 

13.  162 

27.  188 

76.  567 

162.  995 

335.  992 

(7  =  10-2 

0.2 

0.  087 

0.  193 

0.  437 

1.546 

4.  184 

11.347 

.  5 

.761 

2.  555 

6.  735 

17.  823 

1.0 

1.  189 

3.  787 

9.  714 

25.  054 

2.  0 

’.’536 

.943 

1.  907 

5.  701 

14. 121 

35.  229 

.5.0 

1.  232 

1.989 

3.  721 

10.  104 

23.  590 

55.  601 

10.  0 

2.  241 

6.  389 

16.  047 

35.511 

79.  551 

20.  0 

4.  284 

6.  695 

11.  298 

26.  273 

54.  866 

116.  592 

50.  0 

10.  277 

15.  572 

24.  990 

53.  082 

102.  976 

205.  783 

100.  0 

20.  107 

29.  943 

46.  675 

93.  804 

173.  888 

335.  889 

200.  0 

500.  0 

39.  462 

57.839 

87.  843 

168.  225 

300.818 

207.  051 
250.  022 

(7=10-3 

0.2 

0. 169 

0.  309 

0.  613 

1.808 

4.  558 

11.  88.3 

.5 

.  640 

1. 184 

3.  212 

7.  674 

19. 189 

1.0 

’.’729 

1.  151 

2.  024 

5.  103 

11.603 

27.  839 

2.0 

1.  361 

2. 112 

3.  562 

8.  343 

17.  937 

40.  945 

5.0 

3.  256 

4.  913 

7.  886 

16.  749 

33.  .306 

70.  570 

10.  0 

6.  369 

9.  455 

14.  690 

29.  420 

55.  301 

110.  389 

20.0 

12.  5.33 

18.  377 

27.  897 

53. 146 

95.  179 

181.  992 

50.  0 

100.  0 

30.  804 

44.  539 

65.  916 

i  118.842 

202.  064 

9. 181 

200.  0 

65.  632 

9.  896 

500.  0 

--- 

192.  925 

220.  429 

320.  516 

146.  157 

159.  609 

Distance  (kilometers) 

60.6 

1 

1  121 

242 

606 

1,211 

2,  420 
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It  is  believed  that  these  tabulated  values  of  IF 
are  of  considerable  intrinsic  theoretical  interest, 
since  they  can  l)e  used  to  verify  the  validity  of 
operational  methods  used  to  calculate  the  field  of  a 
dipole  near  gently  curved  surfaces.  The  p  series 
formula  for  IF  was  derived  by  such  a  method. 

LTing  eq  (8)  and  (11)  in  conjunction  with  the 
tabulated  values  of  IF,  the  amplitude  E  and  phase 
6  of  the  vertical  electric  field  are  plotted  in  figures 
1  to  16.  The  frequencies  chosen  are  0.2,  0.5,  1.0, 
2,  5,  10,  20,  50,  100,  200,  and  500  kc.  The  ground- 
conductivity  values  are  4,  10~“,  and  10“^  mho/m, 
corresponding  to  sea  water,  well-conducting  land, 
and  poorly  conducting  land,  respectively.  The 
values  shown  in  these  curves  can  be  compared  di¬ 
rectly  with  some  earlier  calculated  results  [1], 
using  essentially  Norton’s  method  [12]  for  the  fre¬ 
quency  range  15  to  500  kc.  The  agreement  is  usu¬ 
ally  within  2  percent  for  the  cases  that  could  be 
compared  directly.  It  appeared  that  the  maxi¬ 
mum  discre])aacies  occurred  in  regions  where  the 
earlier  curves  were  interpolated  between  the  flat- 
earth  values  and  the  first  term  of  the  residue  series. 
As  mentioned  above,  no  interpolation  is  required 
in  the  present  curves,  so  they  are  believed  to  be 
more  accurate. 

To  conform  with  standard  practice  [9],  the  ratio 
k  of  the  effective  to  the  actual  radius  of  the  earth 
is  taken  to  be  4/3  in  the  calculations.  To  illustrate 
the  effect  of  modifying  the  earth’s  radius,  factor  k 
from  1.0  to  4/3  is  illustrated  in  figures  17  and  18. 
The  amplitude  E  and  the  phase  lag  are  plotted 
as  a  function  of  frequency  for  distances  of  124  and 
621  miles  with  a  ground  conductivity  of  10“^  mho/ 
m.  The  difi'erence  between  the  two  sets  of  curves 
for  k=l  and  4/3  is  not  insignificant,  particularly 
from  the  standpoint  of  phase.  However,  at  dis¬ 
tances  less  than  150  miles  or  so,  and  at  frequencies 
less  than  100  kc,  the  amplitude  difference  is  within 
1  percent,  and  the  phase  difference  is  less  than  2 
degrees. 

4.  Conclusion 

It  is  believed  that  the  curves  presented  in  this 
report  will  provide  basic  information  on  ground- 
wave  propagation  in  the  low,  very  low,  and  ultra- 
low  frequency  bands.  It  shoidd  be  em])hasized 
that  the  ionosphere  has  not  been  considered  here, 
and  consequently  the  total  field  at  larger  ranges 


(>200  miles)  can  be  modified.  It  is  possible,  how¬ 
ever,  by  using  special  techniques,  to  eliminate  or 
reduce  the  “sky  wave”  even  at  large  distances  so 
that  the  computed  values  of  the  ground-wave 
propagated  fields  should  have  some  significance 
even  out  to  1,000  miles. 


We  thank  Loris  Peyry  for  assistance  with  the 
calculations,  and  A.  G.  Jean  and  J.  R.  Johler 
for  their  helpful  suggestions. 
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Figure  1. 

Amplitude  of  the  electric  field  as  a  function  of  distance  for  various  frequencies  normalized  so  that  the  radiation  field  is  100  mv/m 
at  Z)  =  l  mile 
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Figure  2. 
See  figure  1  legend. 
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Figure  3. 
S(>e  figure  1  legoiid. 
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See  figure  1  legend. 
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See  figure  1  legend. 
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See  figure  1  Icgcud. 
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Phase  (lag)  of  the  electric  field  as  a  function  of  distance  for  various  frefiuencies  normalized  so  that  the 
phase  is  180°  as  D  approaches  zero. 


See  figure  0  legend. 


Figure  11. 
See  figure  9  legend. 


Figure  12. 
See  figure  9  legend . 
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Figure  13. 
See  figure  9  legend. 


Figure  14. 
See  figure  9  legend. 
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AMPLITUDE  OF  FIELD  STRENGTH  E  (mv/m)  (degrees) 


Figure  15. 
See  figure  9  legend. 


Figure  16. 
See  figure  9  legend 


Figure  17. 
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AMPLITUDE  OF  FIELD  STRENGTH  E  (mv/m) 


Figure  18. 


Boulder,  January  19,  1956. 
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